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I INTRODUCTION 

Physical Vapor Transport (PVT) has been examined numerically for conditions 
where the boundary of the enclosure is subjected to a nonlinear thermal profile. Prior 
numerical studies have (3-9,11 - l3,17 - 20) considered the effects of thermal and solutal con- 
vection, under various thermal boundary conditions, and for conducting (linear tem- 
perature profile) and adiabatic walls. These studies showed that convection is responsi- 
ble for the non-uniformity of the transport flux at the interface. It was also shown that 
for conducting walls, when thermal convection is dominant, the maximum flux occurs 
between 120 ° and 130 0 (depending on the aspect ratio) from the bottom-heated 
configuration. Extremet et al. p3) used a similar set of governing equations to investi- 
gate multizone PVT. They emphasized the importance of solutal convection and 
showed that it can compete with thermal convection to produce multiple cell patterns. 
Zappoli (19) demonstrated that the release of latent heat associated with the evaporation 
and the condensation of the gas at both interfaces can affect the concentration field 
near the interface. However, the actual thermodynamical dependency of the vapor pres- 
sure of nutrient with temperature of the interface is often discarded. The case for mer- 
curous chloride (Hg 2 Cl 2 ) in which the thermodynamical dependency of the vapor pres- 
sure law is taken into account has been studied (26) . It was shown that solutal convec- 
tion is the dominant mode of transport. These studies included heat conduction in the 
crystal, with conducting walls. It was also shown that the interface temperature devi- 
ates from a constant value when a thermal gradient is present in the crystal. Kassemi 
and Duval (24 25) have shown that when internal and surface radiation become important 
in the process, non-uniformities of interfacial mass flux can occur. These past numeri- 
cal studies indicate that under earth gravity, it is difficult to obtain a diffusive flow for 
the horizontal configuration. ^ 

Experimental works in PVT have focused on determining transport rates. Conder 
et al. (1S) have determined experimentally that the transport rate of mercurous iodide is 
controlled by mass transfer in the vapor phase. They also showed that the transport 
rate decreases with the pressure of carrier gas such as argon, and increases with the 


Visiting Research Associate, Hcole Nationale Superieure des Mines, France 



temperature difference between source and sink. In an early work on chemical vapor 
transport (CVT), Chandra and Weidemeier (10) observed the shift of the maximum tran- 
sport rate away from the horizontal configuration, and also showed the transport rate 
increased with the pressure of nutrient (i.e. temperature difference). Related to natural 
convection in PVT, Schiroky and Rosenberger <12) experimentally studied free convec- 
tion of gases in a horizontal cylinder with differentially heated end walls and obtained 
velocity measurements with Laser Doppler Velocimetry. At high Rayleigh numbers, 
pronounced 3-D flow behavior was observed. These measurements were later 
confirmed numerically by Smutek et al. (17) with a 3-D model. 

The goal of this work is to numerically study more extensively solutal convection 
in PVT of mercurous chloride subjected to a nonlinear temperature profile. The effect 
of thermal radiation is neglected. We include the effect of the angle of tilt, and in con- 
trast to the results obtained by Markham and Rosenberger (n) for the case of thermal 
convection, we consider the case when solutal convection is dominant. For practical 
crystal growth applications, it is sometimes necessary to grow the crystal with a non- 
linear thermal profile. We will therefore input a measured thermal profile (27) , with an 
overheat to avoid heterogeneous condensation of the vapor on the walls (as we explain 
in paragraph II). Our studies indicate that this profile destabilizes the flow field when 
thermal convection js dominant for the vertical configuration. However, when solutal 
convection dominates, density gradients associated with solute gradients stabilize the 
flow field for the bottom-heated configuration. Finally, we investigate various cases 
applicable to the microgravity environment, and we deduce the gravity levels required 
to restore a diffusive mode of mass transfer and avoid non-uniformity of the interfacial 
mass flux distribution for the horizontal configuration. 

n APPLICATIONS AND CRYSTAL GROWTH TECHNIQUE 

The mercurous halides Tiave recently become important because of their applica- 
tion in acousto-optic and opto electronic devices. The properties that make this class of 
materials important are: 

a) the large spectral transmission range, 

b) the large photo-elastic coefficient, 

c) the very high acousto-optic figure of merit, " • - ■ 

d) the extremely low acoustic velocity, and 

e) the reasonably low acoustic attenuation. 

These halides belong to the tetragonal (D4h) crystal class and are positively 
uniaxial with large birefringence. The refractive indices for are 2.621 and 1.962 

(birefringent material), the transmission range is 0.36 pm to 20 pm. This wide optical 
transmission range is very unusual and potentially very useful for wideband spectros- 
copic devices. 

The need for high optical quality crystals comes from communication applications 
of the device for which mercurous chloride appears to have great potential. Of the 
metal halides, mercurous chloride has the slowest shear acoustic velocity in the [110] 
direction (347 Also, the acoustic attenuation for high quality crystals appears to 
be reasonably low: 376 dBcm^GHz -2 . These properties are well suited for fabrication of 
Bragg cells which can offer high resolution, high diffraction efficiency, and high 
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dynamic range. 

Mercurous halides decompose to mercury and mercuric halide at the melting 
point ( Hg 2 X 2 -» Hg + HgX 2 ). The existence of the two liquid fields and compositional 
stratification largely preclude the solidification of large Hg^i 2 crystals by conventional 
melt growth techniques. However, mercurous halides exhibit reasonably high vapor 
pressure at low temperature, and crystal growth by vapor transport can be accom- 
plished at rates well within practical range. 

Because of heterogeneous condensation on the walls, a linear thermal gradient is 
sometimes inadequate to grow mercurous chloride crystals. To avoid condensation on 
the walls, practical experimental thermal profiles include an overheat in the vapor 
phase as shown in Figure 1. Unless noted, all simulations reported in this article use 
this thermal profile. 

III. MATHEMATICAL MODEL AND NUMERICAL FORMULATION 

We consider a rectangular enclosure of height H, vapor transport length L, and 
crystal and source material length /, where / = L/2 (Figure 2). The cavity’s axis is 
inclined at an angle 0 with respect to the gravitational acceleration. The solid nutrient 
A exists at both sides of the ampoule in equal amounts. One side is at a temperature 
T s (source), the other at a temperature T c (sink, i.e., the growing interface). With T s > 
T c , the compound A (mercurous chloride) sublimes from the source, is carried with a 
residual gas B (carbon monoxide), and condenses on the crystal interface. The source 
and sink interfaces are assumed to be geometrically flat. Their temperature is deter- 
mined by a heat balance in the gas and the crystals. We impose the temperature profile 
along the upper and lower walls. The thermal conditions at the interfaces impose ther- 
modynamically the weight fraction of component A in equilibrium with the solid 
through the partial pressure of the gaseous phase of compound A: 



We assume complete rejection of B on both interfaces, and bounding walls are 
impermeable to both components. The molecular weights of the two vapor components 
are different, specifically one being mercurous chloride (M A = 472 g/mol), the other 
being carbon monoxide ( M B = 28 g/mol). We selected carbon monoxide because it has 
been found experimentally to be a substantial residual vapor component (2S) . The 
viscosity and thermal conductivity are assumed to be concentration and temperature 
independent. The binary diffusion coefficient is taken to be pressure dependent and 
temperature independent. 

The transport in this enclosure is governed by the system of nonlinear, coupled 
conservation equations for total mass (continuity), balance equation for momentum 
(Navier-Stokes), species mass, and energy. For a two-dimensional enclosure, and 
steady state conditions, the governing equations are the following: 

In the vapor, (/ < x < L+l ; 0 < y < //) ; 

Continuity equation: 
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Conservation of species equation: 
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Conservation of energy in the fluid phase: 
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and in the solid, (0 < x <1 ; L+l < x < L+21 ; 0 < y < //) : 
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Boundary conditions: 

On the upper and lower walls, (/ <x <L+l, y=0 and y=H), the boundary condi- 
tions of impermeability, rigidity and no-slip were used: 


3 a), (x ,0) 3g>, (jc JI) n 

IT s 5 =0 


_v(jc,0) = v(* //) = 0 


(9) 


m(jc,0) = u(xjl) = 0 (10) 

On the source, (x = l ; 0 < y < //), the diffusion of A from the source and impermeabil- 
ity to B were used, in addition to no-slip and an energy balance at the crystal inter- 
face: 
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Analogously on the crystal, (x = L+l ; 0 <><//): 


u(L+l,y) = - 
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The following thermal profile was used as a boundary condition along the ampoule 
(y=0, and y=H), (this equation is expressed in reference to an approximate fit to exper- 
imental data (27) , see Figure 1): 

- 0.10 <; / < 0.22 


563.16 

I T(t) = 608. + 4.97 i- 0.70 l 2 - 5.91 10~ 2 / 3 

7 ~ + 6.67 10"V+ 2.60 10 r 4 / 5 - 2.49 10“5/ 6 

593.16 

At the ampoule ends, the temperatures are prescribed as: 

(*=0 ; 0 <y<H) , T=T S ; (x=2 L ; 0 <y<H ) , T=T C 

Relative to Figure 1, during crystal growth the ampoule (Figure 2) is placed in 
the thermal profile. The hump region corresponds to the location of the vapor com- 
ponent A and B inside the ampoule. The source material lies in the region with the 
larger temperature near r >10cm . Whereas, crystal growth occurs in the region 
corresponding to t<-2cm. In our numerical experiments we position the ampoule in 
the growth region with a temperature less than the source in order to drive the process. 
In addition, the length of the hump region can also be adjusted so that we have a 
much larger source region. With respect to Figure 2, the following transformation is 
used to relate the laboratory reference to the ampoule: 

x = Ki-t 

where Ki is the position of the source and vapor interface in the laboratory reference 
frame. 

An energy balance at the interface between the source material and the fluid 
phase yields equation 13 with a source term accounting for sublimation; likewise a 
sink term (equation 16) also appears for the condensing interface at the crystal end. 
The growth rate of the interface is much smaller than the vapor velocities, typically 
three orders of magnitude. Thus, in comparison to the heat source term and conduction 
terms, the term accounting for latent heat of sublimation is negligible. This is because 
the product term of the growth rate and the latent heat of sublimation is negligibly 
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small. Also, similar to previous investigators we are considering that the interface 
remains flat during the growth of the crystal. 

The system of coupled nonlinear partial differential equations has been solved 
using the SIMPLE algorithm (Semi-Implicit Method for Pressure-Linked Equations 
(30> ). The conservation equations were discretized using a hybrid scheme, which pro- 
vides a second-order spatial accuracy. This two-dimensional version of the algorithm 
uses three staggered grids, one for the scalar variables and one for each of the velocity 
components. In contrast to a vorticity/stream function approach, adopting a primitive 
variable approach allows a straightforward formulation of the velocity boundary condi- 
tions. However, the solution of the continuity equation is no longer automatically 
satisfied. 

For most of the simulations, a 20x60 primary grid was used. The code was run 
on VAX 8800 biprocessor computer, and the iterative process was stopped when a 
Cauchy-Schwartz type criterion was satisfied for each variable. All the results 
presented in this article have been obtained with a convergence criterion of 10"* for 
each of the variables. The code was compared with a more common vorticity/stream 
function algorithm under the same conditions for various Grashof numbers. Results 
provided by the two codes differed by less than one percent (29) . 

IV DISCUSSION 

Residual gas in PVT comes both from the process of sealing the ampoule (carbon 
monoxide) or later degasing of the source material. Because of the inevitable presence 
of this residual gas, or the introduction of a buffer gas, solutal convection can be the 
dominant mode of transport. The effect of the angle of tilt on solutal convection is stu- 
died. We then examine the effect of both partial pressure of carbon monoxide and 
mercurous chloride on the total mass flux. And conclude with a parametric study to 
show the levels of gravity necessary to minimize natural convection. 

AH simulations presented in this part are performed for an ampoule of aspect 
ratio of 10. The dimensionless parameters for the process under the conditions of this 
study are shown in table 1. 


Table 1: Typical dimensionless parameters 

Thermal Grashof number = 

105T 

Solutal Grashof number = 

16500 

Prandtl number = 

0.7 

Schmidt number = 

1 

Aspect ratio = 

10 

1 < Mass Peclet number < 30 | 


IV. 1. Effect of the angle of tilt on the total mass flux in solutal convection 

We now consider the effect of the angle of tilt on the total mass flux when solutal 
convection is dominant. We observe in Figure 3 that, for a solutal Grashof number of 
1650, the maximum transport rate occurs between 30 ° - 50 °. This is in contrast to 
the case for thermal convection at similar Grashof numbers, considered by Markham 
and Rosenberger (!3) where the maximum occurs between 120 ° - 135 ° (depending on 
the aspect ratio), when the cavity is tilted toward the bottom heated configuration. 
There is approximately a 90 ° shift in maximum transport rate between thermally and 
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solutally dominated convection modes. This shift can be attributed to a change in tran- 
sport modes between thermal and solutal convection. Note that the top heated 
configuration is destabilized, since we have heavy fluid over light fluid. Even though 
we use a steady state mathematical formulation, we are fortuitous to obtain a solution, 
only because we are operating below the critical Rayleigh number. For operating 
regions in the neighborhood of the critical Rayleigh number, a time dependent formu- 
lation and solution would be necessary. 

IV.2. Effect of the partial pressure of residual gas on the total mass flux 

For a given thermal boundary condition, which fixes the partial pressure of mer- 
curous chloride at the interfaces through the vapor pressure relationship, and in the 
vertical configuration, we have varied the partial pressure of carbon monoxide at the 
evaporating interface from 0.1 Torr to 100 Torr. The evaporating interface was at a 
temperature of 593 0 K and the condensing interface at 563 ° K. We have taken into 
account the dependence of the diffusion coefficient of carbon monoxide in mercurous 

chloride with pressure (D^ a y). The partial pressure of residual gas affects the tran- 
sport rate because by increasing the partial pressure of carbon monoxide, the mass 
fraction of mercurous chloride at both interfaces decreases. Therefore, the mass flux is 
reduced. Results are shown in Figure 4. The total mass flux decreases exponentially 
with the partial pressure of carbon monoxide. The trend is in agreement with the 
experimental observations of Conder et al. <18) for mercurous iodide (Hgt 2 ) and argon. 

IV.3. Influence of the temperature difference on the total mass flux 

The temperature difference is determined by the location of the interfaces in the 
thermal profile. We selected a transport length of 10 cm between the source and the 
crystal, and an ampoule of aspect ratio 10. A neutral position in the thermal profile 
for which the temperatures of the interfaces are the same can be selected. When the 
ampoule is positioned where the sink temperature is lower than that of the source (see 
Figure 1), the temperature difference drives the transport. Within the covered range 
the temperature difference is proportional to the distance of pulling. Figure 5 shows 
the interface temperature versus the distance of pulling. This figure also shows a sim- 
ple relationship between the temperature difference driving transport and distance. On 
this basis, temperature difference and distance can be used interchangeably. The dis- 
placement of the ampoule from the neutral position is taken as the starting point for 
crystal growth. Thus as the crystal grows, it approaches the neutral point if the 
ampoule is stationary. We have computed the total mass flux versus location of the 
ampoule for three different partial pressures of residual gas. The results are shown in 
Figure 6. We notice that in all three cases the transport rate increases with temperature 

difference following a power law: a &T n with j < n <, j. We also observe 

that, as found in the previous paragraph, the total mass flux decreases with pressure of 
carbon monoxide, resulting in the staggering of the curves. The curves have a steep 
initial slope, where the exponentially increasing concentration difference drives the 
transfer. Then, as shown in Figure 57 as the temperature difference between source and 
sink approaches the neutral position, there is less driving potential for transport due to 
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a decrease in AT, thus the mass flux decreases. These results match the experimental 
results of Conder et al. (18) for a cylindrical ampoule and agree with the 1-D diffusive 
PVT model of Faktor and Garret <2) . 

IV.4. Effect of the temperature boundary condition 

We have previously shown that when a linear axial temperature gradient is 
imposed on the crystal boundary <2fi) , a non-uniform interface temperature distribution 
results. This non-uniform interface temperature can induce recirculating cells in front 
of the interfaces. These results have been obtained for thermal convection with a 
linear gradient of constant sign. However, the experimental thermal profiles in our case 
have gradients of opposite signs at the two interfaces. The presence of the temperature 
peak in the middle of the ampoule destabilizes the flow in both vertical configurations 
for thermal convection. When considering thermal convection, density decreases with 
increasing temperature. Therefore, the density is smaller in the core of the ampoule 
than at the interfaces. As a result, both vertical configurations will be thermally 
unstable. Since our model problem is for steady state conditions, reliable numerical 
results cannot be obtained for these cases because the problem becomes unsteady. 

The introduction of a temperature profile with an overheat throughout the entire 
vapor region does not have much effect on fluid flow pattern in tilted configurations 
for the case when the flow is strongly dominated by solutal convection. When solutal 
convection is dominant, the density gradient resulting from the mass fraction gradient 
stabilizes the flow in the bottom-heated configuration (Figure 7). In this case, the tem- 
perature gradient does not locally reverse the density gradient As a result, the flow 
field remains stable. However, the top heated configuration remains unstable for this 
case. Note again that we were able to obtain a solution for this case only because we 
have a stabilized configuration. 

IV.5. Microgravity environment 

It has been shown that convection results under earth bound conditions for the 
horizontal configuration. Heat conduction in the crystal (26) or internal radiative heat 
transfer (24 - 25) can induce convection even when the gravity vector is parallel to the 
longitudinal thermal gradient. 

To reduce convection, one of the possible alternatives is to grow the crystal in a 
microgravity environment. Configurations with the g-vector perpendicular to the 
ampoule (i.e., horizontal configuration) have been examined for levefs of gravity Tang- 
ing between KT 6 ^ to ig 0 - Figure 8 shows the variation of the total mass flux with 
gravity. We discern two distinct regions. Below and at lO“ 3 go> the flow is diffusive. 
No recirculating cell is present, and the interfacial mass fl ux distribution is jherefore 
planar. Above KT 3 g 0 , the flow switches to convective mode, and a strong recirculating 
cell appears, as shown in Figure 9. The corresponding streamline distribution is in Fig- 
ure 10. At low gravitational acceleration, the Stefan wind drives the flow, and it is 
insensitive to steady levels of gravity less than 10' 3 g 0 - 

V CONCLUSIONS 
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Our study covers the flow field characteristics during the growth of Hg 2 Cl 2 cry- 
stals in a rectangular ampoule under terrestrial and pg conditions for a nonlinear 
prescribed thermal gradient. Taking into account the temperature dependency of the 
equilibrium vapor pressure at the gas/solid interfaces reduces considerably the control 
we have on the process. Solutal convection is the main phenomenon governing the 
fluid flow inside the ampoule. 

The following have been found to be evident: 

1. The magnitude of convection is maximum for the top heated configuration when 
inclined approximtely 40 0 . 

2. The transport rate decreases exponentially with the pressure of residual gas. Pres- 
sures of carbon monoxide varying between 0.1 Torr and 100 Torr have been 
examined. The binary diffusion coefficient was corrected for its pressure depen- 
dency. The numerical trend is in agreement with the experimental results 
observed by Conder et al. (18) . 

3. The transport rate increases with the temperature difference as a power law 
( Jh Si ci 2 a A7" with j <n s y). We have examined three different pressures of 

residual gas. The family of curves indicates a decrease in transport rate as resi- 
dual pressure increases. Conder et al. (18) experimentally found the same depen- 
dency, which is also numerically predicted by the 1-D diffusive PVT model of 
Faktor and Garret (2) . 

4. The experimental temperature profile with an overheat in the core of the ampoule 
applied to avoid heterogeneous condensation of mercurous chloride on the walls, 
apppears to destabilize the flow field when thermal convection is dominant for 
both vertical (0 0 ) top-heated and (180 ° ) bottom-heated configurations. How- 
ever, when solutal convection is dominant, the density gradient resulting from the 
solutal concentration gradient appears to stabilize the flow in the bottom-heated 
configuration. The top-heated configuration becomes unstable for the case con- 
sidered, solutal Grashof number in the neighborhood of 10000. 

5. The microgravity environment provides a means for lowering convection. For 
gravity levels of 10 _3 g 0 or less, the Stefan wind drives the flow. No recirculating 
cell is predicted. 
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NOMENCLATURE 


a 

Ar 


Constant term in the vapor pressure law (a=29.75) 


Aspect ratio. 


L_ 

H 

* 


b Constant term in the vapor pressure law (b=11767.1K) 
C P Thermal heat capacity of the gas 
Dm Binary diffusion coefficient 


F(x,y) source and sink terms. 


dT 


| pc '“V» 


g 

go 

Gr, 


Gr s 

H 

Ja 

J T 

k 

K 

L 

/ 

Mi 

P 

Pe 

Pr 

R 

Sc 

T 

t 

u 

v 


Gravitational acceleration 
Gravitational acceleration at sea level 


Thermal Grashof number. 


Solutal Grashof number. 
Height of the enclosure 




£P 2 // 3 ccAco a 


Mass flux of A 


= p u(L+l,y) 

4 


Total mass flux of A 


J T = fp (L+l ,y)u(L+l ,y)dy 


Thermal conductivity of the gas 
Degree Kelvin 
Vapor transport length 
Length of the crystal 
Molar mass of compound i 
Pressure 


Peclet number, 
Prandtl number. 


In 


PbU'+D 
Pad) l 



Perfect gas law constant ( R = 8.32 kgm 2 s~ 2 mor'K~ i ) 


Schmidt number, 


r ^ 

_ J L .. 

P Dab 


Temperature 

Coordinate attached to the laboratory reference 
Velocity in the x-direction 
Velocity in the y-direction 
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Greek symbols 


a Coefficient of solutal expansion 

P Coefficient of thermal expansion 

k Thermal conductivity of the solid 
p Dynamic viscosity 

Co,- Mass fraction of component i 

p Density of the gas 

0 Angle between the ampoule axis and the gravity vector 
Subscripts 

A Compound A 

B Compound B 

c Pertains to the crystal interface 

s Pertains to the source interface 

T Total 

x Component in x direction 

y Component in y direction 

Superscripts 

c Pertains to the crystal interface 

s Pertains to the source interface 

* Dimensionless 
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Figure 1 .—Temperature profile along the ampoule. 
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Figure 2. — Physical vapor transport enclosure. 
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Figure 3. — Effect of angle of tilt on the total mass flux inside 
the enclosure (Ar = 10, Gr s * 1650, Gr t = 105, top-heated and 
bottom-heated configurations correspond to 0° and 180 c 
respectively). 
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Figure 4.— Effect of the partial pressure of carbon monoxide on 
the total mass flux. (Ar * 10, vertical configuration, Gr s = 1 654 
Gr t = 1050). 
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Figure 5. — Interface temperatures as a function of the distance of 
pulling (L = 1 0 cm, Ar = 1 0). 


Temperature difference (K) 

Figure 6.— Influence of the temperature difference on the total 
mass flux (vertical configuration, Gr s = 16500, Gr t = 1 050, 
Pm - 0.2 Torr, Pm = 1 Torr, P r n ~ 1 0 Torr respectively). 
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Figure 7. — Flow field In vertical configuration 
with experimental thermal profile (Ar = 1 0, 
Gr s m 1 650, Gr t » 1 050). 
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Figure 8. — Total mass flux as a function of the gravity level 
(horizontal configuration, 1 O^gg < g < g^. 



Figure 9.— Flow field for various gravitational levels (g 0 , 10** 1 g 0 , 
1 0 _2 g 0 , lO^gQ correspond to Gr s = 16500, 1650, 165, 16.5 
respectively). 
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Figure 10. — Corresponding streamlines (g 0 , KHqo.IO^q, IO-^q 
correspond to Gr* = 16500, 1650, 165, 16.5 respectively). 
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